function [ AResp, PhResp, Omega, w ] = GetExactPhiOneElPlate( w, R_p, a, Bi)
%[AResp, PhResp, Omega, w] = GetExactPhiOneElPlate(w, R_p, a, Bi)
%   
    Omega = zeros(size(w));
    AResp = zeros(size(w));
    PhResp = zeros(size(w));
    for i=1:length(w)
        Omega(i) = w(i)*R_p^2/a;
        x=sqrt(Omega(i)/2);
        a1=sin(x);
        a2=sinh(x);
        a3=cos(x);
        a4=cosh(x);
        a5=a3*a4 + (a2*a3-a1*a4)*x/Bi;
        a6=a1*a2 + (a1*a4+a2*a3)*x/Bi;
        a7=a1*a2*a3 + a3*a4*a5;
        a8=a1*a2*a5 - a3*a4*a6;
        AResp(i)=20*log10(sqrt(a7^2+a8^2)/(a5^2+a6^2));
        PhResp(i)=atan( a8/a7 );
   end
    
end

